AD-A227  330 


Naval  Research  Laboratory 


Washington,  OC  20375-5000 


NRL  Memorandum  Report  6616 


me  hLc 


^uh1/ 


The  Approach  to  Equilibrium  in  a  Reflex  Triode 


Robert  e.  Terry  and  John  M.  Les 

Radiation  Hydrodynamics  Branch 
Plasma  Physics  Division 


James  L.  Geary 
Berkelev  Research  Associates 


September  24,  1990 


DTIC 

ELECTE I 

OCT  1 2 1990 1 


This  research  was  sponsored  by  the  Defense  Nuclear  Agency  under  Subtask  Code  and  Title:  RL 
RA/Advanced  Technology  Development,  Work  Unit  Code  and  Title,  00079,  Advanced  Technology 
Development.  This  work  was  supported  under  MIPR  Number  89-565. 


Approved  for  public  release,  distribution  unlimited 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OM8  No  0704  0188 


Public  rfOOnmg  burden  tor  thi\  collection  Of  information  'S  CMimated  TO  average  1  *Our  oef  '*\oO''se  including  the  lime  <0'  reviewing  instructions  searcning  exiting  data  SO«cc« 
gathering  and  maintaining  the  data  needed,  and  competing  and  reviewing  me  collection  of  information  Send  comments  regarding  this  Durden  estimate  or  any  other  asoect  of  this 
collection  of  information,  including  suggestions  »or  reducing  this  burden  10  Washington  neadauarte's  Services.  Directorate  for  information  Operations  and  Reoons.  >2  IS  teMerson 
Oavis  Highway.  Suite  1204.  Arlington,  va  22202-4302  and  to  the  Office  o<  Management  and  Budget.  Paperwork  Reduction  Project  (0  704-0  ’88)  Washington.  DC  20503 


1.  AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  OATE 

1990  September  24 


4.  TITLE  ANO  SUBTITLE 

The  Approach  to  Equilibrium  in  a  Reflex 
Triode 


3.  REPORT  TYPE  ANO  OATES  COVERED 

Interim 


FUNDING  NUMBERS 

PE  -  62715H 
WU  -  DN880-191 


6.  AUTHOR(S) 

Robert  E.  Terry,  John  M.  Les 
and  James  L.  Geary* 


7.  PERFORMING  ORGANIZATION  NAME(S)  ANO  AOORESS(ES) 

Naval  Research  Laboratory 
Washington,  DC  20375-5000 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


NRL  Memorandum 
Report  6616 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  ANO  AOORESS(ES) 


10.  SPONSORING  /MONITORING 
AGENCY  REPORT  NUMBER 


Defense  Nuclear  Agency 
Alexandria,  Virginia  22310 


11.  SUPPLEMENTARY  NOTES 


♦Berkeley  Research  Associates 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution 
unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

^  Analytic  and  PIC  descriptions  are  combined  to  provide  a  theory  of  reflex  triode  operation  in  and 
near  the  steady  state.  The  spectrum  of  reflecting  electrons  is  examined  both  for  a  steady  input  of 
electrons  into  the  scattering  foil  and  for  the  dynamically  evolving  process  of  achieving  the  reflexing 
equilibrium.  Simple  approximations  to  analytic  modifications  of  the  classical  theory  then  offer  a  more 
versatile  theory  of  the  re  flexing  system.  Kr  .  .•  7'*'1-  «  •  '  ,f  Lt:c'’r,~ 


14.  SUBJECT  TERMS 

Reflex  triode 
Vlasov  theory 


17.  security  classification 
OF  REPORT 

UNCLASSIFIED 


NSN  7540-01-280-5500 


Electron  scattering 
PIC  simulation 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


15.  NUMBER  OF  PAGES 

33 


16  PRICE  CODE 


20.  LIMITATION  OF  ABSTRACT 


standard  form  298  (Rev  2-89) 

Pi^c'ioro  by  ANSI  5td  2  39  <8 
298  >02 


CONTENTS 


INTRODUCTION  . 

PHENOMENOLOGY  OF  THE  REFLEX  TRIODE  . 

BASIC  RELATIONS  FROM  THE  VLASOV  PICTURE 

STEADY  STATE  REFLEX  TRIODE  MODEL  . 

AUXILLARY  MODELS  . 

SIMULATING  THE  APPROACH  TO  STEADY  STATE 

CIRCUIT  MODELS  FOR  TRIODE  DEVICES . 

CONCLUSIONS  . 

ACKNOWLEDGEMENTS  . 

REFERENCES  . 

DISTRIBUTION  LIST  . 


1 

3 

4 
6 
12 
16 
18 

22 

23 

25 


THE  APPROACH  TO  EQUILIBRIUM  IN  A  REFLEX  TRIODE 


Introduction 

Unique  in  the  spectrum  of  opening  switch  concepts,  the  reflex  triode1,2  offers 
the  possibility  of  triggering  the  opening  phase  to  control  the  overall  synchronism  of 
a  inductive  store  pulser.  The  basic  triode  construction  (and  even  its  more  elabo¬ 
rate  variations )  has  an  attractive  simplicity  when  compared  to  the  complex  source 
hardware  for  plasma  erosion  opening  switches,  density  controlled  opening  switches 
or  plasma  filled  diodes.  In  common  with  the  plasma  filled  diode,  the  reflex  tri¬ 
ode  switch  is  its  own  load  and  thus  minimizes  the  inductance  ” downstream ”  of  the 
opening  switch,  a  property  long  known  to  be  helpful  in  any  inductive  store  concept. 

There  are  two  principal  concerns  in  bringing  such  switch  into  routine  use.  A 
relatively  minor  problem  is  the  requirement  for  an  axial  magnetic  guide  field,  which 
leads  to  auxiliary  coil  and  shielding  hardware.  The  more  critical  problem  is  the  fact 
that  the  larger  currents  are  not  presently  achievable,  nor  credibly  implied  by  the 
heretofore  successful  scaling  laws  concerning  the  ’’anomalous  gap  closure ”,  which 
limits  the  peak  conduction  phase  current. 

The  physics  that  determines  the  upper  current  limit  for  the  conduction  phase 
is  not  fully  understood  at  present.  The  theory  presented  here  seeks  (i.)  to  establish 
a  more  complete  physical  basis  for  reflex  triode  switch  operation,  (ii.)  to  refine,  in 
the  context  of  real  devices,  some  of  the  ciassicai  notions  held  from  early  equilibrium 
theory  and,  (Hi)  to  provide  some  conceptually  simple  models  which  allow  estimates 
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particle  approach  readily  incorporates  sufficient  generality  to  treat  the  global  switch 
behaviour  and  to  provide  qualitative  guidance  as  to  the  relative  importance  of  var¬ 
ious  phenomena.  The  continuum  theory  approach  readily  captures ,  quantitatively, 
(i)  the  observed  current  density  enhancements  in  the  conduction  phase  and  (ii)  the 
typically  observed  opening  phase  current  and  voltage  levels  -  even  though  we  use 
very  simple  ideas  to  model  the  opening  process. 

A  diagram  of  a  typical  reflex  triode  is  shown  in  Figure  1,  annotated  with  the 
usual  variables  for  the  familiar  Child-Langmuir  diode  theory  and  some  new  ones. 

REFLEX  TRIODE  SYSTEM  VARIABLES 

=  o  ,  9'  =  0  9  =  1  ,  9'  =  y  9  =  90  ,  9'  =  9'0 

Ki  fi  K2 


Here  Th  is  the  reflexing  electron  flux  (unidirectional),  ritTl  are  the  ion  and 
neutral  fluxes  generated  at  the  anode,  ¥(£)  =  V(£)/V0  and  'P'(£)  =  ( ^ ) ^~  arc  the 
normalized  potential  and  field  profiles.  The  boundary  conditions  at  each  electrode 
are  indicated.  The  unsubscripted  7  factor  is  a  measure  of  the  mean  field  at  the 
anode  surface  required  to  extract  an  ion  -  due  to  charge  exchange  drag.  Drag 
factors  in  the  ’’vacuum”  region,  also  discussed  below,  are  subscripted  7  variables. 
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for  the  effects  of  ” anomalous  gap  closure”  or  alternate  hypothetical  opening  mech¬ 
anisms  on  the  predicted  device  waveforms. 

This  requires  that  we  address,  using  a  quasi-steady  state  Vlasov  formulation, 
and  a  dynamic  PIC  formulation,  some  important  questions  which  arise  in  connection 
with  the  reflex  triode  switch: 

•  Does  charge  exchange  with  neutrals  provide  the  rapid  gap  closures  observed; 
how  can  one  derive  and  improve  the  observed  scaling  laws? 

•  What  are  the  fundamental  limits  on  the  current  density ?  Are  they  determined 
by  the  electrode  physics,  microinstabilities,  or  electron  neutral  collisions? 

We  will  first  build  an  extension  of  the  ” classical ”  theory  to  develop  those  ideas 
needed  in  connecting  the  steady  state,  single  particle  Vlasov  metaequilibria  to  slow 
timescale  changes  in  the  boundary  conditions  and  reflexing  spectra.  We  will  illus¬ 
trate  some  connections  and  contrasts  between  the  PIC  treatment  and  the  continuum 
theory.  Next  we  will  show  the  development  of  a  circuit  model  which  addresses  the 
initial  buildup  phase  of  the  current  and  the  transition  from  low  to  high  impedance. 
Using  some  simplified  results  of  the  more  complicated  models  we  will  compare  the 
circuit  model  picture  to  experimental  data. 

Phenomenology  of  the  Reflex  Triode 

Present  reflex  triode  switch  devices  are  generally  operated  on  the  early  phase 
of  the  current  pulse  discharging  energy  into  the  storage  inductor.  The  triode  first 
presents  a  high  voltage  to  the  driver  as  the  electron  population  builds  up  and 
begins  to  oscillate  ( reflex )  through  the  foil  anode.  Those  electrons  transiting  the 
primary  gap  (the  so-called  A-Kl  gap)  at  high  energy  are  stopped  by  the  secondary 
cathode  and  begin  to  charge  it  negative  relative  to  the  anode.  This  process  continues 
until  there  is  a  detailed  balance  between  received  and  re-emitted  electron  current 
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and  a  floating  potential  is  established  on  the  secondary  cathode  K2.  This  floating 
potential  usually  exhibits  one  or  more  positive  voltage  steps  during  the  conduction 
phase,  probably  due  to  bursts  of  ion  current  from  the  anode.  The  ion  current  is,  of 
course,  unidirectional  -  no  ions  are  re-emitted.  Yet,  so  long  as  there  exists  a  voltage 
difference  on  the  A-K2  gap  which  is  large  enough  to  boost  an  electron  re-emitted 
from  K2  to  a  one  range  energy  in  the  foil  anode,  viz.  >  10 kV ,  the  reflex  electrons 
can  continue  to  populate  a  spectrum  of  energies  above  that  required  for  capture 
and  thereby  enhance  the  ion  emission  from  the  anode. 

It  is  the  ion  emission  that  is  the  key  to  the  low  impedance  state  which  then 
develops  after  a  few  tenths  of  ps.  For,  while  the  ion  current  is  relatively  small  (par¬ 
ticularly  for  the  true  triode  geometries),  the  ion  space  charge  in  the  gap  neutralizes 
the  primary  electron  space  charge  and  shortens  the  effective  gap  spacings  near  each 
electrode. 

The  low  impedance  state  would  persist  forever  (in  principle )  but  always  falls 
prey  (in  practice)  to  some  catastrophe  associated  with  the  movement  of  the  floating 
conductor  nearer  the  anode  potential.  The  opening  phase  occurs  too  soon  in  the 
driver  pulse  to  offer  effective  power  transfer  and  fails  to  follow  previously  observed 
conduction  time  scalings,  based  on  the  notion  that  charge  exchange  neutrals  fill  the 
triode  and  short  the  gaps. 

Basic  Relations  from  the  Vlasov  Piet ure 

In  contrast  to  the  highly  discrete  particle  code  approach,  discussed  later,  the 
treatment  presented  first  relies  upon  the  continuous  representation  of  the  reflex  tri¬ 
ode  plasma  by  a  single  particle  distribution  function.  Both  techniques  have  serious 
limitations  and  the  relative  merits  of  either  can  be  discussed  at  length,  but  the  best 
use  of  such  diversity  appears  to  lie  in  cross  comparison.  Briefly  stated,  the  discrete 


The  electron  energy  spectrum  is  g(a)  (1  >  a  >  'P,,)  and  arises  through  multiple 
scattering  events  in  the  anode  foil.  In  these  terms  two  basic  variables  are  the  bipolar 
Child-Langmuir  current,  I0  =  2.334  •  10~®  j|-Vo3^2(l  +  \/ Zme/m,i),  and  the  ratio  of 
ion  to  electron  current,  jij  je  —  a v/Zmf/mt. 

Electrons  scattered  in  the  anode  with  each  reflexing  passage  populate  a  dis¬ 
tribution  fe  which  can  be  projected  onto  the  subspace  [W0  <g>  p)  of  (gyrokinetic) 
adiabatic  invariants3 

W-pBz  =  l-mtVz2-e<I>  =  W0  =  -eV0a  , 
mew2  c  \\ 

m  =  0  -  c(v*,*))  • 

These  variables  form  a  "good”  set  because,  even  though  ji  is  only  approximate,  it 
will  be  ’’used”  only  for  one  electron  transit  in  and  out  of  the  random  scattering 
process,  since  we  fiduciate  the  constants  of  motion  to  the  anode  plane.  The  capture 
energy  at  the  anode  is  thus  mapped  to  the  limit  W0  — >  -~eV0  (cr  =  1),  while  the 
limit  W„  — >  0  (<r  =  0)  confers  the  full  A-Kl  voltage  to  an  primary  beam  electron. 
The  modification  to  ji,  due  to  the  largely  axial  acceleration  v  is  obtained  only  in 
the  limit  that  the  magnetic  field  due  to  the  conduction  phase  current  is  competitive 
with  the  guide  field  in  the  switch. 

The  ’’free  streaming”  electron  populations  then  admit  a  projection  operation 
II w„fe  —  /e(x,v |<r)  such  that 

je( x)  =  eT(i)  =  e  J  =  e  J  d(T  [  J  d3v  w,/e(x,  v|<r)  ]  , 

separating  the  reflexing  current  into  subcomponents.  The  subcomponents  may  be 
considered  either  as  discrete  6a  or  differential  da  sections  of  the  distribution  function 
defined  on  subsets  of  the  [x,  vj  domain.  The  important  property  is  that  the  Vlasov 


5 


equation  is  trivially  satisfied  for  any  subcomponent  fe(x,v\a),  in  the  free  streaming 
region. 

A  formal  derivative  with  respect  to  W0  or  a  defines  the  usual  spectrum  of 
free  electron  energies,  =  J  d3vYlw0fe,  which  must  satisfy  two  a  priori 

constraints. 


(i.)  The  integral  over  all  energies  must  provide  a  measure  of  the  total  number  of 
reflexing  electrons.  With  Te  =  Jv  >Q  dzv  vzfe(z  =  0,  v|tr  =  0)  ,  the  emitted  source 
flux  at  Kl,  the  normalization  used  here  is 


rRefux(<T)  =  f°  da !  =  vre  G(a)  :G(  1)  =  1  . 

Jo  deTi 

(ii.)  The  gradient  with  respect  to  the  energy  ratio  a  must  represent  a  positive 
definite  electron  probability  density 


dr 


Re  flex 


da 


dG  f 

=  :  g{<r)  =  —  =  /  d3vvzUWofe  >0 

dcT  Jv,>  0 


Any  time  dependence  of  the  reflexing  process  or  the  injection  flux  Tc  is  reflected 
mostly  in  the  time  development  of  TfTe  and  to  a  lesser  extent  in  the  spectral  form 
factor  g.  It  will  be  shown  that  very  small  changes  in  tj  near  some  critical  value  can 
make  enormous  changes  in  Te. 

Steady  State  Reflex  Triode  Model 

In  order  to  make  contact  with  some  of  the  phenomena  seen  in  the  reflex  triode 
switch  one  must  extend  the  steady  state  theory 4  of  Prono,  Creedon,  et.al.,  to  include 

i.  more  general  emission  boundary  conditions,  e.g.  charge  exchange  presheath 
effects; 

ii.  asymmetric  potential  distributions  and  gap  spacings; 

Hi.  realistic  foil  scattering  spectra  (of  parallel  energy); 
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iv.  electron  and  ion  drag  mechanisms; 

v.  deterministic  opening  mechanisms,  viz.  changes  in  the  internal  metaequi¬ 
librium  state  which  raise  the  impedance  quickly. 

Only  this  sort  of  extension  can  provide  a  simple  test  bed  for  experimental  compar¬ 
ison  that  incorporates  the  results  of  the  PIC  code  REFLEX  whenever  possible. 

The  usual  method5  for  examining  space  charge  limited  flow  is  to  determine  the 
velocity  of  a  current  carrier  as  a  function  of  the  local  potential  value.  Adopting 
the  simple  notion  that  the  ” average ”  charge  carrier  may  be  subject  to  some  modest 
drag  forces  -  whether  from  collective  effects  or  direct  collisions  -  a  ” test  particle” 
variation  on  this  theme  can  be  constructed. 


Let  T(£)  represent  the  kine4'c  energy  of  a  particle  relative  to  its  free  accel¬ 
eration  value  over  the  full  gap,  viz.  mv2/2qV0.  When  the  drag  force  is  modeled 
as  proportional  to  the  kinetic  energy  divided  by  a  ” mean  free  path  for  momentum 
transfer”  or  ’’damping  length”  (X),  the  equation  of  motion  can  be  integrated  once 
to  yield 


dT  ^ 

li  ±  2lT  = 


where  7  =  d/X  and  the  prime  also  denotes  d/d(.  The  damping  length  scale  is  now 
just  a  parameter  that  can  be  obtained  from  auxiliary  models.  Denoting  the  various 
drag  processes  by  jiie  ,  the  kinetic  energy  becomes 


Ui) 


-2 


+ 


e27'  *'(6) 


> 


and 

Ti{()  =  +  J  e27i  U-6)  *'(^)  . 

In  the  limit  that  V  •  J  =  0  ,  the  kinetic  energy  terms  can  represent  quasi-steady 
fluxes  of  electron  and  ion  space  charge  in  the  Poisson  equation.  The  first  integral 
can  be  obtained  to  provide  an  integral  equation  formulation  for  ^'(0 
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*'(« 


=  [  y.F(*,*„)  +n,/2«)  +  + 

7.  f  i?.7;1/3({.)  -  «7i  f(  l'7' 

Jo  Jo 


where  cold  emission  has  been  assumed  and  the  spectral  contribution  T  is  given  by 


F(9,9o) 


‘ 


The  integral  over  the  spectrum  can  be  simplified  -  via  Liebniz  rule,  a  partial  inte¬ 
gration  and  a  reversal  of  integration  order  -  to  give 

^(*,¥0)  =  g0(*-*o)3'  ~  /  <%'(*)(¥ -0)’  . 

j'fo 

The  constraint  ’P  >  is  understood;  T  vanishes  elsewhere. 


In  the  limit  of  zero  damping,  strict  space  charge  limited  ion  emission,  and 
symmetric  potentials  this  formulation  reduces  to  that  of  Creedon.  If  furthermore, 
the  reflexing  flux  is  set  to  zero,  there  obtains  a  full  reduction  to  bipolar  Child- 
Langmuir  theory. 


When  one  can  simplify  this  integral  equation  because  the  modifications  due  to 
the  drag  effects  are  very  small,  viz.  7^  1  ,  then  the  integrated  terms  above  reduce 

to  simple  mean  vaiu  ?  estimates.  On  the  [A,K1]  domain  the  normalized  electric  field 
is  then  approximately 


=  fa  kv^  +  a(exp-7„/2)[(V(l-<W-l)]  +  ¥„)]  * 

O  d  j  </(*) 


Here  the  effects  of  damping  mechanisms  are  approximated  by  de  and  d 


de  =  \/(l  -  27e(l  ~  exp  — 7e)  , 


d{  -  (exp  7i)(l  -  27i)  +  , 

where  the  7 ’s  represent,  as  defined  above,  the  reciprocal  mean  (dimensionless)  free 
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path  for  moment  1  n  transfer  and  the  reductions  to  the  classical  formulation  are 
preserved.  To  stay  sensible,  these  rough  forms  involving  the  damping  factors  must 
be  confined  to  the  domains  7e  <  0.864  and  7j  <  0.5 

The  continuity  of  determines  the  normal  K2  electric  field  parameter,  which 
figures  in  the  gap  breakdown  dynamics,  'P'0('I'0,£0,a,/?),  as 

*'o  =  ldeV*~o  +  a(r7,J-,de,(ii,7)(v/(l-£ft^0)-l)]i  . 

Usually  the  parameter  tj,  which  measures  the  overall  strength  of  the  reflexing 
flux  relative  to  the  primary  flux  from  the  cathode  (K1 ),  is  considered  to  be  a  function 
of  the  voltage  on  the  primary  gap  and  thus  increases  in  rough  proportion  to  the 
number  of  reflexing  passes  this  voltage  would  imply  for  a  given  foil  anode.  In  an 
actual  triode,  allowing  for  the  re-emission  of  electrons  from  the  floating  conductor, 
this  simplification  no  longer  holds.  In  general  the  re-emitted  population  is  bounded 
above  in  energy  by  the  value  of  the  floating  potential  difference  1  —  'J'„.  As  -*  1, 
particularly  in  the  case  of  pure  space  charge  limited  ion  emission  (7  —*  0),  the 
eiectron  population  is  held  to  such  low  energies  that  the  flux  ratio  rj  can  never  be 
much  larger  than  1  even  if  electrons  reflex  10  times  or  more  through  the  anode. 
If  one  defines  <7r(l,'I,0)  =  /(l,  ,P0)^/1  —  then  the  largest  sdlowed  value  of  tj  is 
constrained  by 


Vmax  _ 


<  - 


-Si 


-A 


9/1672 

4/  1  1  —  \/(l  -  V^)(l  +  y^)  '  v'l-* 


+ 


in  order  to  retain  a  real  solution  at  the  entry  point  to  the  ( spatially )  bounded 
reflexing  region  in  the  primary  gap. 

Only  because  we  have  approximated  the  integral  equation  character  of  the 
problem  can  the  A-Kl  current  be  determined  by  the  familiar  quadrature 


fi  =  3  f1  d9 

-Wo  g\(0) 


69 


ln(4g+/go) 

qT 


Here  0  measures  the  current  density  enhancement  over  a  Child-Langmuir  bipolar 
diode  -  due  to  the  reflexing  population  as  measured  by  the  dux  ratio  q.  In  exactly 
the  same  way  as  the  ” classical ”  formulation,  this  quadrature  implies  essentially 
singular  current  in  the  allowed  metaequilibrium  states  if  the  parameter  q  approaches 
a  critical  value.  Such  a  critical  value  exists  for  all  <  1,  viz.  low  impedance  states 
persist  until  the  absolute  value  of  the  A-K2  voltage  can’t  get  an  electron  to  a  1  range 
energy  in  the  anode.  Even  though  qmax  can  drop  close  to  1  for  — »  1  there  always 
exists  a  value  qcrit  <  qmax  which  offers  the  singular  state.  The  logarithmic  nature 
of  this  singularity  has  been  indicated  in  the  formula  above,  where  Q0  represents  the 
true  minimum  value  of  the  (dimensionless)  electric  field  profile  function  and  G+  it’s 
value  on  the  boundaries  of  a  domain  66  in  the  potential.  Under  what  conditions 
real  experiments  can  approach  this  sort  of  equilibrium  is  simply  unknown,  but 
the  laboratory  experience  does  indicate  that  a  high  level  (0  >  1000 )  of  current 
enhancement  is  obtained  routinely.  From  the  theory  presented  here  this  observation 
implies  that  the  interior  electric  field  is  being  shielded  to  better  than  a  few  parts 
in  107  relative  to  its  peak  value  near  the  electrode  surfaces.  Indeed  the  modest  0 
field  structure  is  no  different  from  that  characterizing  very  large  0  values.  Figure 
2  shows  a  'Po  =  0.5,  0  =  224  triode  field  profile.  In  this  case  the  value  of  q  differs 
only  slightly  from  the  critical  value  at  which  0  — *  oo.  Those  few  signifcant  figures 
in  q  correspond  to  charge  variations  of  at  least  several  hundred  million  electrons 
per  cc  in  the  denser  space  charge  regions,  well  above  thermal  fluctuation  levels. 

The  simulation  experience  confirms  the  global  robustness  of  this  low  impedance 
state,  but  does  not  resolve  the  large  current  density  enhancements  one  observes. 
The  particle  model  is  limited  by  random  noise  in  its  ability  to  suppress  the  interior 
field  values  on  its  grid.  Moreover,  the  several  hundred  million  electrons  variation 
noted  above  corresponds  to  only  a  Jem  code  macroparticles.  In  resonance  with  the 
equilibrium  picture,  the  particle  code  does  not  exhibit  0  values  much  in  excess  of 
200! 
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Reflex  Trinde  Field  Profiles 


On  the  other  hand ,  the  genera/  character  of  the  simulation ’s  Vlasov  phase  space 
is  in  good  agreement  with  a  similar  picture,  Figure  3,  drawn  from  the  equilibrium 
theory,  showing  the  two  beam  populations  and  the  contours  of  constant  cr,  again  for 
the  same  parameters  as  in  the  previous  figure.  The  probability  density  /e(x,  v |cr)  is 
also  constant  on  each  of  these  contours  if  there  is  a  true  steady  state,  but  in  any  case 
these  contours  are  the  Vlaxov  characteristics  for  the  problem.  The  primary  feature 
to  note  here  is  that,  for  a  real  triode  with  'P „  >  0,  the  hole  in  phase  space  between 
the  beam  and  reflexing  populations  will  be  depleted  by  the  circuit  constraints  and 
filled  by  turbulence  in  the  primary  beam  and  reflexing  electron  population.  When 
and  at  what  level  the  balance  is  struck  will  be  shown  in  the  PIC  simulations. 

It  may  be  argued  that  such  a  delicate  balance  (±108  electrons)  is  not  physically 
accessible  to  any  real  system  and,  moreover,  that  inferences  drawn  concerning  the 
operating  point  and  the  sensitivity  of  this  model  to  boundary  condition  changes 
could  not  therefore  be  expected  to  hold  under  laboratory  conditions.  Ideas  of  this 
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Reflex  Triode  Phase  Space 

contours  of  constant  electron  phase  density,  g(<?) 


nature  also  require  proof,  however  appealing  they  may  be  intuitively.  In  this  regard 
one  fact  must  be  kept  in  mind  -  the  sensitivity  of  the  ” nearly  singular”  equilibria  is 
logically  inseparable  from  their  ability  to  explain  the  large  current  density  enhance¬ 
ments.  If  these  sensitivities  and  other  properties  are  unsatisfactory,  then  another 
explanation  for  the  enhanced  current  density  must  be  proposed  and  examined.  The 
only  model  that  explains  the  current  density  enhancement  is,  at  present,  the  ’’clas¬ 
sical”  one  or  the  improved  version  discussed  above.  Because  the  equilibrium  model 
exhibits  sensitivity  to  boundary  conditions,  interior  Held  extinction  levels,  or  dissi¬ 
pative  processes,  then  experiments  can  and  should  be  designed  to  test  these  ideas 
insofar  as  it  is  practical  to  do  so. 

Auxiliary  Models 

The  Vlasov  descriptions  just  discussed  require  as  input  several  nontrivial  bits  of 
physics  from  either  experiment  or  other  theoretical  treatments.  Chief  among  these 
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is  a  scattering  model  for  the  effective  spectrum  of  electron  energies.  The  next  most 
critical  is  an  estimate  for  the  flux  of  fast  charge  exchange  neutrals  Tn  produced  at 
the  anode  and  probably  responsible  for  the  rapid  gap  closure  phenomena.  Finally 
one  must  consider  the  magnitude  of  various  drag  mechanisms  on  both  ions  and 
electrons,  the  local  field  strengths  and  neutral  pressures  near  K2  which  will  figure 
in  the  gap  shorting  time,  and  the  boundary  condition  evolution  at  the  anode. 

Aside  from  the  steady  conduction  phase  evolution  oft]  to  values  near  the  sin¬ 
gular  point,  one  finds  that  (from  an  operational  point  of  view  in  fixing  (3)  all  such 
auxiliary  physics  changes  ( time  variable  anode  boundary  conditions ,  reflexing  pop¬ 
ulation  changes,  or  damping  processes )  are  imbedded  -  via  the  spectrum  T  and  the 
various  j’s  -  in  the  single  parameter  a(Tj,iF,de,di,y) 

(1  +  i./2)i.  +  ff(  1  -  *„)>/»  -  9/167! 

°  e-»/*[i_(i_7j/2)(l -*)■/»] 

which  determines  the  relative  level  of  ion  current  in  the  triode  by  forcing  the  bound¬ 
ary  condition  on  the  electric  field  at  the  anode. 

Generally  as  a  increases  from  1  the  ion  space  charge  brings  0  up,  and  conversely. 
Among  the  factors  that  lower  0  by  reducing  a  are  (i)  any  reduction  in  ion  emission 
caused  by  anode  resistance  7,  ( ii )  any  movement  of  towards  1,  and  (Hi)  any 
reduction  in  the  parameter  tj.  The  direct  effect  of  the  approximate  damping  factors 
on  0(a)  is  to  increase  it,  presumably  due  to  the  longer  time  spent  by  a  particle  in 
the  gap.  On  the  other  hand  one  expects  the  exponential  character  of  the  damping  to 
come  into  play  as  the  7 terms  exceed  their  present  constraint  values  and  lower  a. 
In  addition,  the  drag  of  slow,  reflexing  electrons  on  slow  moving  ions  and  neutrals 
late  in  the  conduction  phase  would  be  manifested  in  a  direct  lowering  of  q  by  pushing 
more  electrons  to  low  energy  domains.  A  good  evaluation  of  the  drag  effects  requires 
the  solution  of  the  full  integral  equation. 
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Evaluating  the  spectral  convolution 


One  way  to  obtain  a  spectrum  is  to  use  REFLEX’s  Moliere 6  scattering  model 
which  operates  in  the  limit  where  the  mean  squared  scattering  angle  is  larger  than 
the  usual  Rutherford  result.  Here  the  (relativistic)  electrons  exhibit  a  toted  scat¬ 
tering  cross  section  of  the  general  form:  <tt  oc  Z1^3(Z  -f  l)/i/2.  In  such  a  mylar 
foil,  for  example,  a  lOkeV  electron  has  a  2  micron  range.  The  scattering  model  has 
been  benchmarked  with  the  TIGER  code,  c.f.  J.M.  Les  and  J.  Ambrosiano1 . 

The  spectrum  of  energies  g(o-)  can  be  determined  over  the  full  reflexing  voltage 
range,  and  used  to  set  T .  In  Figure  4  an  overlay  of  many  such  flux  spectra  shows 
an  essentially  universal  shape  for  all  the  voltages  examined  ( V0  —*  lOOfcF ),  only  the 
normalization  changes  and  the  tj  parameter  holds  this  information  in  any  case.  A 
power  law  fit  (a  ss  0.68)  is  observed 

*  (i 

and  this  result  allows  an  analytical  evaluation  of  the  spectral  convolution.  One 
finds,  with  p  =  |  ~ , 

W».)  =  (1  +  a)VT^  [,»/»  +  y,»/’,r1(,-*>,|/>)] . 

This  spectrum  is  the  asymptotic  result  for  the  pure  scattering  model  where  the 
input  flux  (at  <r  =  0)  is  held  fixed  and  no  collective  plasma  effects  are  active  in 
changing  the  distribution.  Modifications  of  this  spectrum  are  discussed  below  in 
connection  with  the  PIC  code  simulations. 
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Figure  4 

Evaluating  the  ion  charge  exchange  sources 

A  Monte  Carlo  code 8  has  been  developed  which  treats  the  ion  charge  exchange 
process  through  a  given  background  of  neutral  gas;  the  benchmarks  and  early  appli¬ 
cations  are  encouraging  -  the  correct  spatial  dependence  of  neutral  flux  is  predicted 
for  the  case  of  constant  cross  section.  The  input  required  is  a  neutral  density  spatial 
profile,  an  electrostatic  potential  function,  or  a  flux  of  hot  ions,  and  an  appropriate 
cross  section.  The  output  is  a  neutral  flux  profile  and  energy  distribuion  which  can 
be  used  as  a  first  step  in  calculating  the  scaling  laws  for  the  gap  closure.  Some 
simple  estimates  for  the  velocity  of  charge  exchange  neutrals  have  been  obtained 9 . 
While  these  are  of  limited  generality,  the  code  and  analysis  would  predict  final  neu¬ 
tral  velocities  in  the  range  of  20  — *  60cm/ps  for  the  expected  field  strengths  in  a 
typical  /3  %  1000  system.  This  corresponds  to  the  acceleration  of  hydrogen  ions 
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over  a  potential  of  a  few  hundred  volts  to  perhaps  a  kilovolt  in  a  small  fraction  of  a 
centimeter.  We  also  find  a  weakening  dependence  of  this  neutral  velocity  with  ion 
layer  thickness.  After  the  layer  is  thickness  (xD)  is  large  enough  to  get  ions  up  to  the 
rollover  energy  of  the  charge  exchange  cross  section  (W0  ~  lOfceV )  in  the  ambient 
electric  field,  viz.  when  qEx0  %  Wot  then  the  final  velocity  is  very  weakly  depen¬ 
dent  on  the  ambient  E  field.  Moreover,  for  the  thinner  layers,  the  neutral  energy 
scales  linearly  with  the  electric  field  -  in  agreement  with  rough  analytic  scaling  laws 
and  the  empirical  scaling  rule  tj  J3  oc  d^K\  •  At  these  velocities  the  observed  time 
difference  t$  between  the  onset  of  opening  and  the  ’’anomalous  closure”  of  the  main 
gap  -  with  the  consequent  prompt  lowering  of  the  triode  voltage  -  is  accounted  for 
by  the  difference  in  gap  spacing.  For  typical  bank  and  gap  parameters  one  would 
expect  this  time  difference  to  be  about  0.10  1 1. 

Simulating  the  Approach  to  Steady  State 

Using  the  REFLEX  code  one  can  determine  the  spectrum  and  flux  ratios  that 
arise  in  a  triode  system  as  it  comes  into  a  fully  developed  reflexing  state.  As  the 
electron  population  settles  down  and  establishes  the  two  sheaths,  the  spectrum 
relaxes  as  well,  showing  a  marked  smoothness  in  the  time  averaged  population  at 
any  particular  value  for  the  <r  parameter.  The  time  average  is  carried  out  over  ten 
or  more  electron  transit  periods,  as  estimated  from  the  local  simulation  parameters. 

As  a  direct  comparison  with  the  Vlasov  description,  the  secondary  circuit  is 
isolated  resistively  so  that  a  fixed  fraction  of  the  primary  gap  driver  voltage  tends  to 
build  up  over  time.  The  steady  state  achieved  tends  to  oscillate  about  a  secondary 
voltage  value  greater  than  the  primary,  but  collective  effects  can  force  rather  wide 
excursions  about  this  mean. 

If  the  secondary  voltage  is  greater  than  the  primary  (\P0  <  0),  then  the  re 
emitted  electrons  from  K2  can  easily  flow  back  into  the  primary  gap  although  most 
of  the  population  is  turned  around  without  ever  striking  K2.  A  reflexing  state  is 
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achieved  with  a  rather  broad  phase  space  measure  in  vz  near  the  anode,  and  the  time 
averaged  spectrum  is  peaked  near  a  =  0  and  is  monotonically  decaying  elsewhere. 

For  (Vo  >  Oj  the  reflexing  population  is  narrower  in  vz  near  the  anode  and 
the  quality  of  reflexing  is  improved  in  that  more  primary  current  is  drawn  and  the 
voltage  is  further  reduced  on  the  primary  gap.  As  shown  in  Figure  5a  where  contours 
of  equal  particle  population  are  plotted  on  the  time  and  energy  domain,  the  time 
history  of  the  energy  constant  is  dominated  by  the  detailed  wave  particle  interaction. 
The  phase  space  hole  is  energetically  possible  but  gets  filled  in  rather  completely 
by  microturbulence.  The  potential  profile  (Figure  5b)  is  shown  for  K2  well  above 
the  Kl  potential,  but  the  phase  space  is  dearly  populated  in  the  region  where  the 
gap  should  be,  as  is  evident  in  Figure  5c,  where  the  time  averaged  spectrum  g(cr )  is 
shown.  The  spectrum  shows  a  more  pronounced  plateau  near  cr  ^  but  the  region 
near  <r  ~  0  is  heavily  populated  by  particles  of  essentially  full  primary  gap  energy. 
The  scattering  model  may  be  responsible  here  insofar  as  the  Moliere  limit  is  going 
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Figure  5a 
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to  be  progressively  violated  as  the  electrons  lose  energy  -  a  more  robust  population 
for  moderate  cr  values  would  probably  increase  the  Child-Langmuir  enhancement. 


POTBNTLAL  AT  T  =  D.2660E-06 


SPECTRUM  G  RT  T=2 .  5R90E-O7 


Figure  5c 

Circuit  Models  for  Triode  Devices 

One  can  apply  a  circuit  model  to  this  triode  theory  when  uj  <C  wpe,  and  it  is 
worth  noting  some  basic  estimates  for  the  electron  density  and  plasma  frequency  in 
the  switch.  The  following  are  lower  bounds  appropriate  near  the  anode  where  the 
electrons  move  most  rapidly: 
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-  -  W-3'10'")^)^)2. 


These  timescaJes  for  the  response  of  electrons  to  changes  in  the  ion  charge  distribu¬ 
tion  and  the  electrode  surface  charge  make  a  very  reasonable  case  for  a  quasi-static 
picture.  While  the  timescedes  observed  in  REFLEX  to  recover  a  low  impedance 
state  from  some  disruption  can  refine  this  constraint  further ,  the  robustness  of  the 
low  impedance  states  observed  in  the  particle  code  argues  for  the  viability  of  the 
near  equilibrium  picture.  The  a  priori  stability  of  the  system  could  not  be  assessed 
from  any  equilibrium  model,  but  the  basic  result  is  that  the  global  field  and  density 
profiles  demanded  by  the  equilibrium  constraint  V-J  =  0  are  in  fact  preferred  states 
of  the  time  evolving  system  as  well.  The  issue  of  microstability  is  quite  a  different 
story  and  here  both  Vlasov  and  particle  models  speak  volumes  as  to  the  turbulent 
channels  available  to  the  system. 

Whatever  the  details,  if  one  makes  the  assumption  of  near  equilibrium  states, 
then  the  value  of  j3  determines  the  triode  impedance  and  requires  a  \'ak\  large 
enough  to  carry  the  current  supplied  by  the  generator,  viz. 

i'  ^akx  a,  wl.a 

V  SWITCH  =  77 -  =  • 

’oanfc  P 


This  implies  a  single  ODE  for  T,  which  is  of  the  form 

fr  “ 

with  t  =  t/ti  and  1  =  I /hank)  and  =  —  T  .  A (r)  is  a  measure  of  the  initial 
operating  (reflexing)  voltage  and  is  time  dependent  only  because  the  ’’anomalous 
gap  closure"  is  presumed. 


19 


V switch  is  normally  steady  In  the  conduction  phase,  any  process  in  the  triode 
which  forces  a  lower  (3  value  precipitates  the  opening.  In  modeling  the  opening  phase 
it  is  important  to  increase  the  speed  of  the  ion  motion  as  the  inductive  voltage  is 
applied.  One  finds  too  slow  a  decay  in  the  current  enhancement  if  the  ions  just  drift 
at  the  speed  given  by  the  reflexing  voltage.  If  /?*  is  the  peak  current  enhancement 
just  prior  to  opening,  then  the  widening  of  the  gap  between  the  anode  and  the 
neutral  plasma  mass  will  force  the  effective  (3  down  according  to  the  square  of  the 
gap  ratio.  The  application  of  this  notion  in  the  context  of  our  circuit  relation  results 
in  a  time  dependent  f3  given  by 


where  u,  %  21.5  is  a  measure  of  the  ion  drift  speed  and  £0  is  a  measure  of  the  initial 
gap  fraction  covered  by  neutral  plasma.  Folding  all  the  results  together  one  can  get 
a  fair  picture  of  the  conduction  phase  and  the  opening  phase,  induced  by  the  quick 
drop  of  f3  as  the  ions  are  swept  out  of  the  gap.  The  time  delay  for  the  onset  of 
opening  is  determined  experimentally,  only  more  detailed  theory  can  pin  down  such 
a  number.  The  initial  bank  charge  is  calculated  from  experimental  data,  as  well 
as  the  initial  value  of  (3  at  the  observed  onset  of  reflexing.  After  that  the  model  is 
on  it  own,  equipped  with  an  estimate  for  the  anomalous  gap  closure  time  from  the 
charge  exchange  calculations  mentioned  above.  The  peak  current  obtained  and  the 
peak  voltage  at  opening  are  not  pre-ordained  in  any  way  except  by  the  uniqueness 
of  the  ODE  solution. 

The  application  of  this  circuit  model  to  a  typical  shot 10  yields  the  waveforms 
shown  in  Figure  6.  Here  the  peak  current  predicted  is  54  kA  while  the  experiment 
saw  56  kA;  the  predicted  peak  voltage  is  634  kV  while  the  experiment  saw  560  kV. 


20 


Model  Mavcfcrws:  HOSES  I  shot  2887 


Conclusions 

When  viewed  in  the  context  of  the  REFLEX  particle  code  work  documented 
elsewhere  in  these  proceedings,  the  qualitative  picture  of  the  reflex  switch  which 
arises  is  both  simple  and  compelling.  The  device  operates  at  an  enhanced  current 
density  because  neutral  plasma  plays  the  role  of  a  conducting  medium  that  effec¬ 
tively  shortens  the  initial  gap  spacing  to  two  smaller  gaps  (for  a  triode)  as  reflexing 
electron  flux  builds  up.  Each  basic  notion  of  the  ’’classical”  Reflex  Switch  model 
can  be  partially  substantiated  by  these  diverse  theoretical  approaches  and  by  exper¬ 
iment.  Taken  as  a  whole  these  results  show  the  device  to  be  based  on  solid  physics, 
but  clearly  not  understood  well  enough  (as  yet )  to  be  certain  it  can  fullflll  its  in¬ 
tended  role.  We  have  detailed  three  points  of  sensitivity  possessed  by  the  reflexing 
state:  ion  charge  exchange,  electron  drag  from  either  turbulence  or  neutral  interac¬ 
tion,  and  ion  emission  boundary  conditions.  Each  of  these  may  operate  along  with 
the  ” classical ”  opening  mechanism,  viz.  the  shorting  of  K2  to  A,  in  determining 
the  opening  time  and,  hence,  the  maximum  current.  Several  modifications  to  the 
” classical ”  picture  are  implied  by  our  new  formulation,  in  particular 
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•  Vcrit  is  coupled  to  and  no  longer  is  purely  a  function  of  voltage; 

•  High  current,  low  impedance  equilibria  persist  until  the  A-K2  potential  is  less 
than  a  range  energy  «  lOkeV ; 

•  The  triode  phase  space  is  especially  unstable,  but  the  K2  potential  will  work 
to  keep  the  hole  open; 

•  The  spectral  convolutions  are  simplified  and  can  often  be  performed  analyti¬ 
cally,  yielding  quantitative  agreement  with  the  observed  current  enhancement 

(3- 

We  have  found  that  the  reflexing  state  is  potentially  susceptible  to  the  starva¬ 
tion  of  ion  current  due  to  charge  exchange  or  the  decay  of  reflexing  electron  flux  due 
to  electron  neutral  collisions.  If  such  sensitivity  does  indeed  occur  in  the  labora¬ 
tory,  this  would  result  in  z  rapid  increase  in  impedance.  No  experiments  yet  fielded 
can  preclude  any  of  these  ideas  but  these  mechanisms  can  be  quantified  with  some 
effort  and  can  help  to  understand  the  opening  process.  Moreover,  these  additions 
can  be  tested  in  a  relevant  and  cheap  circuit  model  to  gain  insight  and  to  suggest 
future  experiments.  The  circuit  model  clearly  operates  around  the  experimentally 
observed  current  and  voltage  traces  -  without  benefit  of  any  free  parameters. 
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